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We present a channel model useful for analyzing the effects of 
multipath fading in digital radio systems. The frequency response of 
a fading channel is represented by the function A — uBi + juAi, 
where u {=2mf) is measured from the center of the channel, and A , 
A\, and B\ are variable coefficients that change slowly with time. The 
model consists of this function, the joint probability density function 
(pdf) for the three coefficients, and the average number of seconds 
per heavy-fading month for which this response applies. The model 
is derived from a large base of multipath fading data, obtained on a 
26.4-mile path in Georgia in June 1977. It consists of nearly 25,000 
recorded measurements of received power vs frequency in a 25.3-MHz 
bandwidth at 6 GHz. In this paper, we present the methods of data 
reduction and statistical analysis used to derive the model; describe 
some assessments of its validity; and discuss its limitations, virtues 
and possible uses. By all available measures, the model is highly 
accurate. It suffers from a potentially important phase ambiguity that 
can be resolved only via new, coherent measurements. The existing 
model should prove very useful in the design and planning of such 
measurements. 

I. INTRODUCTION 

Multipath fading in microwave digital radio systems can be a major 
source of outage and, therefore, has been the subject of numerous 
recent investigations. 1 " 4 One important objective of current activity is 
to develop a statistical model of fading useful for estimating its effect 
on specific systems and for indicating possible methods of correction. 

Two multipath fading characterizations have been proposed recently 
for terrestrial systems in the channelized common carrier bands. One 
represents the multipath fading frequency response as that due to an 
equivalent three-path medium. 5 The result is a three-parameter, com- 
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plex gain function, where the joint probability distribution for the 
three parameters describes the fading statistically. The second char- 
acterization, which is the subject of this paper, expresses the multipath 
fading frequency response as a complex polynomial expanded about 
the channel center frequency. 6 Specifically, if 

complex channel gain at any time ,.v 

= complex channel gain during nonfading ' 

where a(=2irf) is measured from the center frequency, then the 
polynomial representation is 

N 

H c (io) = A + 2 (A„ +jB n )(ja) n , (2) 

n=l 

where the A- and S-coefficients vary slowly relative to the speeds of 
typical digital radio systems. 

Note that H c {0) is just A , a real number, i.e., we arbitrarily (and 
with no loss in generality) specify the phase shift at the channel center 
frequency to be zero. 

The polynomial function of eq. (2) can fit frequency responses of 
any shape with arbitrary accuracy merely by choosing N sufficiently 
large. Moreover, the difficulty of fitting coherently measured responses 
does not grow significantly with N since H c (u) is linearly related to the 
characterizing parameters (the A n 's and B n 's). 

Several features make the polynomial representation attractive for 
multipath modeling. One is that it leads to simple methods for analyz- 
ing signal processing (note that (ja) a corresponds to the nth time 
derivative). A second is that it suggests a promising idea for adaptive 
equalization; specifically, the rational function 1/H c (w) may be easy to 
realize adaptively when the complex zeros of H c (u), (2), have negative 
real parts. A third feature is parsimony, suggested by the conjecture 6 
that a first-order function (N = 1) may be sufficiently accurate for 
terrestrial paths in the channelized bands below 15 GHz. In that case, 
the channel response function could be characterized by three linear 
coefficients (A , Ai, and Bi). 

The work reported here has confirmed this conjecture for the data 
base used by Rummler. 5 The central activity has been to compute A , 
A\ , and B\ from the recorded measurements and to analyze their joint 
statistics. The outcome is a model for the joint pdf, p(Ao, A\ , B\ ), that 
can be used to assess multipath effects in any digital radio system. 

Section III presents a mathematical description of the new model. 
Section II gives an overview of the data reductions and analyses 
leading to this model, while Section IV gives an expanded discussion 
for readers interested in more details. 

Section V discusses the accuracy of the derived model and presents 
supporting evidence for the steps taken to simplify it. Part of this 
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evidence is obtained by calculating outage times for numerous system 
approaches, using model versions of differing complexity. This exercise 
not only validates the simplest version but illustrates its important 
application to system studies. Section VI summarizes the uses and 
limitations of the model and cites the need for further measurements 
to reinforce and improve it. 

II. DATA REDUCTIONS 

2. 1 The data base 

The multipath fading data base is a set of measurements made on 
a 26.4-mile path in Palmetto, Georgia in June 1977. 5 During fading 
events, received power was measured at each of 24 frequencies distrib- 
uted about 6034.2 MHz and spaced by 1.1 MHz. The measuring filter 
at each frequency had a 200-KHz bandwidth. The power samples were 
normalized by the unfaded received power level (as estimated from 
mid-day measurements), converted to decibels, quantized in 1-dB 
steps, and recorded on magnetic tape. Because of an equipment failure, 
data for the 19th frequency in the sequence (7.15 MHz above the 
center of the channel) were not obtained. 

Each tape record therefore consists of 23 power response values, in 
quantized decibels, corresponding to a single measurement. No phase 
response data were recorded, resulting in unresolvable ambiguities in 
the reductions, as discussed later. Measurements were recorded at 
rates of either five per second or one per two seconds, depending on 
the speed of change in the fading response. The data base used here 
consists of 24,920 records spanning, roughly, 8400 seconds of fading. 

2.2 The polynomial coefficients 

The main reductions consist of (i) determining the coefficients in 
eq. (2) for each record in the data base and (ii) finding a mathematical 
description of the joint probability density function (pdf) of these 
coefficients. In this subsection, we outline the strategy used for deter- 
mining the coefficients. 

Each data record consists of the quantized values of —10 logio | H c \ 2 
at 23 frequencies, where H c is the complex channel gain. 

We denote the recorded values by the set {Pi}, where i is the 
frequency index (i = 1, 24 excluding 19). The data reduction per record 
begins by converting each decibel quantity P, into a power ratio 
Pi(— 10~ p,/1(l ), and fitting the sequence of p, vs frequency with an Mth- 
order polynomial, 

q(u) - Do + uDi + • • • w m D m + • • • w M D M . (3) 

In obtaining the coefficients (Do, • • • Dm), a form of least-squares 
optimization has been used, and M-values of 0, 2, and 4 have been 
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tried. We have evaluated the errors between the p,-'s and the fitting 
function q(u) for each record and have studied the statistics of these 
errors over the ensemble of records that comprise the data base. Our 
finding is that, given the noise statistics and bandwidth of the mea- 
surements, the most suitable polynomial order is M = 2. For, when M 
= 0, the fitting function lacks the curvature needed to model the true 
fading response, while, when M = 4, the third- and fourth-order 
polynomial terms serve more to fit the measurement noise than to fit 
the underlying response. The curvature provided by q(u) when M = 2 
seems to be just right for the 25.3-MHz measurement bandwidth. For 
that case, we estimate the rms fractional difference between q(u) and 
the true power response to be 5 percent at all frequencies. 

Given a value for M and a procedure for obtaining the D m 's from the 
data records, we must choose the order (N) of the polynomial in eq. 
(2) and derive the A n 's and B„s in that equation from the D m 's. Before 
deciding on N, we note that a power-gain function can be obtained 
from eq. (2) having the form 

| H c (a>) | 2 = Do + coA + • • • u> m D m + • • • u 2N D 2N (4) 

where the An's are simply related to the A n 's and B n 's. The objective 
is to match this function to q(u) in eq." (3). 

One possibility is to choose N = Af/2 and to find the A„'s and B n 's 
for which D — Do, A = 2>i, etc. The problem is that, if D M happens 
to be negative, this important term cannot be matched using all real 
An's and B„'s. The approach we have used instead is to specify N = M, 
and to choose the (2AT + 1) A n 's and B n 's so that 

D ={ D ™> °^ m ^ N (5) 

Um \ 0; N<m<2N. w 

This procedure, however, is limited in at least one and sometimes two 
respects. The fundamental limitation concerns those A„'s and J3„'s 
contained within the imaginary part of H c {u), i.e., (A\, A 3 • • • £2, B 4 
• • • )• Given a valid solution to eq. (5), reversing all these coefficients in 
sign results in exactly the same function for | H c (w) | z . Thus, there are 
two possible solutions, each of which is equally valid given the data. 

The polarity ambiguity in the subset (Ai, As, ••• Ba, B 4 •••) is 
reflected as a polarity ambiguity in the phase response associated with 
i/ c (oj), an inevitable consequence of noncoherent measurements. Our 
reduction approach is to assume equally likely polarities for any given 
fade. In computing A n 's and J3„'s for each record, we have randomly 
selected the polarity in accordance with this assumption. 

A second limitation on the application of eq. (5) arises when the 
Dm's for n > N cannot all be forced to zero using real A n 's and B n 's. In 
performing reductions for N = 2, for example, a solution with real A n 's 
and fi„'s is only possible when 
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A = (D 2 - D?/4Do) > 0. 



(6) 



We will limit all remaining discussions to the case N = M = 2. 

The condition A > corresponds roughly to q(u) being concave over 
the measurement bandwidth. In such cases, A 2 = B 2 = 0, and the 
function for H c (io) reduces to a first-order polynomial. When A < 
{q(w) more or less convex over the measurement bandwidth), no real 
combination of A„'s and 5„'s will satisfy all of eq. (5). This condition 
is found to occur over 42 percent of the data base. Our strategy for 
these cases is to choose (Ao, A u A 2 , Bi, B 2 ) so that (Do, Di, D 2 ) 
matches (Do, D u D 2 ); D 3 = 0; and D 4 (= Ai + Bl) has the smallest 
magnitude possible for real A„'s and B n 's. The result is a fourth-order 
component in \H c (i>>) | 2 , namely, u 4 D 4 , which we have found to be 
relatively small (in a sense that we define later). 

The above strategy leads to the following set of solutions: 



A = y/D ; 

S, = -D 1 /2Ao; 

±VA; A>0 



A,- 



A 2 = 



0; A > 

Af/- . Ai 

Ao 



A<0 



^O + gi I; ^<« 



and 



(7a) 
(7b) 



(7c) 



(7d) 



(7e) 



B 2 = A 2 B 1 /A i . 

In computing these coefficients, the polarity of A i is chosen randomly 
for each record. 

The question arises of how the model should take account of the 
two distinct conditions on A. One idea is to treat all computed com- 
binations of (Ao, Ai, Bi) as part of a single statistical population, 
without differentiating between the conditions A > and A < 0, and to 
assume that A 2 = B 2 = in all cases. The result would be a first-order 
polynomial for H c (u) defined by a single joint pdf of Ao, Ai and B x . 
This approach has been adopted for reasons of modeling simplicity, 
and is justified by data given in Section V. 

2.3 The coefficient statistics 

We begin the statistical modeling with the coefficient A , represent- 
ing the complex channel gain at to = 0. The probability distribution of 
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Ao is found to be approximately log-normal, i.e., if 

ADB A 20 logioAo, (8) 

then the pdf of ADB is approximately Gaussian, with a mean (/a) of 
—21.39 dB and a standard deviation (2) of 6.562 dB. For convenience, 
we shall deal with ADB in terms of the standardized variable 

a A(ADB-(i)/-Z. (9) 

The precise pdf for a a used in the model is presented in Section III. 

The joint probability law for Ao, Ai, and B\ can be represented by 
the joint pdf of A\ and B\, conditioned on a a , times the pdf of a Q . 
Furthermore, system study results (Section V) indicate that A\ and 
Bi, for given a Q , can be modeled as independent variables. Hence, we 
write 

p(Ai, Bi | a Q ) - Pa{Ai I ao)pi)(Bi \ a ). (10) 

To model pa( | ) and pa( | ), the entire population of data records was 
divided into 11 subpopulations, each corresponding to a specific inter- 
val of a . For example, all records with a between —2.25 and —1.75 
constituted one such subpopulation, all those with a between —1.75 
and —1.25 constituted another, and so on. Within each subpopulation, 
we computed the mean and standard deviation of A\ , as well as its 
probability distribution, and similarly for B\* 

Both A\ and Bi are found to be essentially Gaussian in every 
subpopulation; in every case, moreover, each variable has a zero mean 
and a standard deviation that varies with the mid- value of a . Conse- 
quently, the joint pdf of A\ and B\, conditioned on a , reduces to a 
product of two Gaussian pdf 's, each with a standard deviation that is 
a function of a Q . By finding suitable mathematical expressions for 
these functions, the statistical modeling of Ao, A\ , and B\ is completed. 

A final quantity to determine is the number of seconds, Tm, for 
which the multipath fading response applies in a heavy-fading month. 
The data base used here represents about 8400 seconds of fading 
which, Lundgren and Rummler have estimated, corresponds to two- 
thirds of a heavy-fading month for the measured path and frequency. 4 
In attempting to relate this finding to other paths and frequencies, we 
have assumed that Tm is proportional to the multipath occurrence 
factor developed by Barnett. 7 This assumption permits the general 
expression for Tm given in the following section. 



* The choice of 11 subpopulations and the particular a -intervals used is somewhat 
arbitrary. Our aim was to obtain both good resolution (narrow intervals of Oo ) and high 
estimation accuracy (many samples per subpopulation). 
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III. THE MODEL 

In this description, all numerical quantities in braces, { } , are data- 
derived constants that may vary with path length, antennas, locale, 
year, etc. It is hoped that the underlying model structure, consisting of 
the form of the transfer function and the functions for the coefficient 
pdf 's, is generally applicable and that the constants alone might be 
subject to change. In any event, what follows is an estimated fading 
model based on reductions of the given data. 

(i) The complex transfer function of a channel, normalized by its 
unfaded gain, is 

1 1 + y'O during non-fading periods; 

H c (u) = A — uBi + juAi during T M seconds per heavy- 
fading month. (11) 

(ii) By assuming T M to be proportional to the multipath occurrence 
factor, 7 we obtain 

T M = {0.11} cFd\ (12) 

where c is the terrain factor, F is the system frequency in gigahertz, 
and d is the path length in miles. 

(Hi) The joint pdf of A , Ai , and Bi (where A is dimensionless and 
A\ and B i are in units of seconds) can be represented by 

p(a Q , Ai, Bi) =p A (A x | cio)pb(Bi | a )p a (a ), (13) 

where 

20logioA - {-21.39} „., 

°° = {6^62} * U4) 



{iv) The pdf of a a is the nearly Gaussian function 



p a (a ) =—=exp 

V27T 



--[do + z{a )f 



dz(a„) 



(15) 



da,, 

where z(a ) is a small nonlinear term in a a given by 

z(a ) = {0.0742}ai; + {0.0125}ai (16) 

If this component were zero, p a (a ) would be precisely Gaussian. 
Figure la shows p a (a ) and compares it with the Gaussian pdf having 
the same mean (0.0) and variance (1.0). The difference seems small 
but has been found to be significant, as we will discuss in Section 5.2. 
(v) The conditional pdf's of Ai and B u p A {Ai \ a a ) and p B (Bi \ a a ), 
are both Gaussian with zero means and with standard deviations given 
by 

a A = Max{{0.14}, [{0.309} + {0.13} a a ]} X 10 -9 s (17) 
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Fig. 1 — Functions used in the multipath fading model: (a) the pdf of a„, compared 
with a normal (Gaussian) pdf; and (b) the standard deviations of Ai and B\, conditioned 
on a,,. 

and 

a B = Min{{0.24}, Max({0.12}, [{0.18} + {0.046}a o ])} X 10" 9 s. (18) 

These functions, which model the variations of a a and ob with a , are 
shown graphically in Fig. lb. They are derived by estimating oa and 
ob from finite numbers of records in nine a -intervals of finite width. 
The results are therefore quite approximate; our tests indicate, how- 
ever, that they are not sensitive to the precise positions or widths of 
the a -intervals. 

IV. DATA REDUCTIONS: EXPANDED DISCUSSION 

This section presents mathematical, numerical, and graphical details 
underlying the data reductions and the construction of the model. The 
interested reader can thereby scrutinize the various stages of analysis 
and reasoning that led to the results in Section III. The less interested 
reader can, with no loss in continuity, proceed to Section V. 
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4. 1 Data organization and notation used 

The data base consisted originally of 24,920 records of decibel power 
gain vs frequency. An early decision was made to remove the first 100 
records in the sequence, as they contained a number of normal (non- 
fading) responses. The model is derived from the resulting set of 24,820 
records. 

The second decision made was to create four subsets of data records, 
each containing about 2160 records distributed randomly over the 
total set (24,820 records). The aim was to verify the statistical regular- 
ity of the data base, i.e., to establish that no handful of records 
dominates the statistics and that individual subsets are representative 
of the data base as a whole. Having verified this regularity (by 
comparing results among the four subsets and between the four subsets 
and the total set), we were able to perform certain costly reductions 
using a relatively small set of records. 

All computations in this study were performed using, as a data base, 
either the four individual subsets, the four subsets merged together 
(8640 records), or the total set of 24,820 records. In every case, we have 
found it helpful to regard the prevailing data base as an ensemble of 
records, each record being a sequence of 23 power gains. The notations, 
symbols, and definitions used in all computations are summarized in 
Table I, which should be consulted throughout subsequent discussions. 

4.2 Deriving the polynomial coefficients 

The process begins with 24,820 records of P, vs i and ends with a 
population of coefficient sets, (A , Ai,Bi), whose reduction to a joint 
probability law (Sections III and 4.3) defines the fading model. Three 
stages of computation are involved here, namely, (i) adjusting the 
recorded P,'s to remove systematic calibration differences across the 
23 measurement frequencies, («') determining, from these adjusted 
data, the optimal polynomial order and the resulting record-by-record 
coefficients of the power gain function, eq. (3), and (Hi) determining, 
from these coefficients, the record-by-record coefficients of the com- 
plex gain function, eq. (2). We now discuss these three stages in turn. 

4.2. 1 Data adjustments 

From physical considerations, we would expect the ensemble prob- 
ability law for Pi to be the same for all i. Thus, the power gain at any 
given frequency should exhibit the same statistical behavior as that at 
any other frequency. This expectation is confirmed by the data, which 
show that the probability law for Pi is approximately Gaussian for all 
i. At the same time, however, the data show small but discernible 
variations of P, and o, with i. 
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Table I — Symbols and Definitions 



Symbol Definition 



i Index of frequency ( i = 1 for lowest frequency; i = 24 for highest fre- 

quency). 

Xi Average of x, over ensemble of records comprising data base ( x, can be any 

variable). 

( Xi ) Average of x, over i for given record; 

1 24 

i?, Actual channel power gain, in dB, at ith frequency. 

P, Recorded channel power gain, in dB, at ith frequency. 

p, Power gain ratio at ith frequency;/?, =10 p,/l °. 

q(u) Polynomial fitted to sequence of p,'s; q(u) = £ D m u m . 

m-0 

Qi Decibel value of q(u) at ith frequency (represents fitted approximation to 

Pi); Qi = 10 log 9 (»,-) ■ 10 log qi. 

of Ensemble variance of Pf, ol =~P1 - (ft) 2 . 

e, Measurement error, (Pi — i?,), at ith frequency; e, = oLm, all i. 

Et Observed error, (P, - Qi), at rth frequency. 

Ci Modeling error, (Qi — i?,), at ith frequency. 

E'ina Mean-square observed error per record; E^ms =(Ej). 

E 'r mB Approximation to £? ms (quadratic in qi); 

_ A 10 pi-qj 

£rms " \llnl0 P i 

Ki, Bias in observed errors per record; Eb = {Et) 

E\ k) Ath root of Ath central moment of E t over ensemble; E[ k) = V(Ei - £,)*. 

These variations were computed for the total record set and for each 
of the four record subsets. Results are given in Fig. 2, where the 
vertical bars span values for the four subsets and the solid curves 
connect the values for the total set. Also given are results for the 
decibel value of p,. The variations for this quantity are virtually 
parallel to those for P, , a consequence of the fact that the pdf of P, has 
roughly the same shape and variance for every i. 

For all three quantities (P,, a, and pi), the consistency of results is 
evident among the four subsets and between those subsets and the 
total set. For each quantity, moreover, the average variation with i is 
more pronounced than the spread among the data subsets. These 
findings make clear that the variations with i are not statistical but, 
rather, the result of systematic calibration differences in the measure- 
ments. The variations with i of P,, in particular, can be explained in 
terms of such systematic differences. Comparable explanations for the 
variations of a, are not as forthcoming. Nevertheless, a decision was 
made to attribute the nonuniformities of both P, and o, to systematic 
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effects, and to compensate for them using simple adjustments of all 
the Pi's in the data base. The adjustments are of the form 

Pi^diPi + Pr, eachi, (19) 

where the variations of a, and /?, with i are chosen to serve two 
purposes. One is to render the ensemble means and variances of the 
adjusted P,'s uniform over i. The other is to conserve the '"global" 
mean and variance (P and a 2 ) of Pi, i.e., the mean and variance of the 
Pi's taken over the entire data base (24,820 records X 23 samples/ 
record). It is easy to show that this is achieved by choosing a, and fit 
as 

«,=-, 0i = P--Pr, eachi'. (20) 

Oi Oi 

The variations of P if a,, and 10 logp, for the entire record set, when 
the data are adjusted according to (19) and (20), are given in Fig. 2 by 
the dashed curves. We see that, by making P, and a, uniform over i, 
10 log p, has been made nearly uniform also, as we expect it to be from 
physical considerations. The a,'s and /?,'s that produce these results 
are close to 1.0 and 0.0 dB, respectively, for all i* 

Though minor, the data adjustments have measurably reduced the 
curve-fitting errors associated with the next stage of reduction. This is 
made evident later. 

4.2.2 Coefficients of the power gain polynomial 

4.2.2.1 Least Squares Method. We wish to compute, for given M, 
those D m 's in eq. (3) that yield the best approximation to the channel 
power response in each record. Ideally, we would like to mmimize 
some measure of (Q, — Ri) (see Table I), such as its mean-square 
average over i, (e 2 ) . Since the values of P, are unknowable, the next 
best thing is to minimize some measure of (P, — Qi), such as (Z? 2 ) = 
E^ms. However, the Q,'s are logarithms of samples of q(u), and so a 
precise least-squares minimization is not analytically tractable. For- 
tunately, so long as | Ei \ < 1 dB for most i, £L can be accurately 
approximated by i? 2 ™*, as given in Table I. This quantity lends itself to 
least-squares minimization. 

For a given M and a given record, the following approach is used to 
find Do, • • • Dm: (i) the 23 P,'s are converted to power ratios, /?,; («) 
the function E^os is expressed in terms of the D m 's in eq. (3); (Hi) the 
M + 1 derivatives of ^^ with respect to Do, • • • D M are computed 
and equated to zero; and (iv) the resulting M + 1 linear equations in 
Do, • - - Dm are solved using matrix methods. 



* The a,'s range from 0.945 to 1.029, and the fr'a range from -0.99 dB to 1.25 dB. 
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As a refinement to this procedure, we add one more step: The 
derived D m 's are multiplied by a common factor chosen to force the 
bias error, (E,) = Eb, to zero.* Thus, for example, if the initial solution 
for the D m 's leads to a bias error E B = X, then adjusting each D m by a 
multiplying factor 10* /10 will reduce Eb to zero. This step reduces E rms 
and makes the Q's unbiased estimators of the P,'s. 

Figure 3 shows two sample records for P„ taken from the data base. 
In each case, the function Q(u) = 10 log q{u), as computed using the 
above procedure, is plotted for M = 0, 2, and 4. Note that, for M = 0, 
the fitting function lacks the curvature needed to accommodate typical 
fading patterns. For M = 4, on the other hand, the fitting function 
provides added curvature that may be artificial, i.e., it may serve to fit 
the measurement noise more than the underlying response. We now 
turn to a discussion of fitting errors and the most suitable choice for 
M. 

4.2.2.2 Random Error Model. To begin, assume that 23 noisy samples 
of Pi are to be fitted by Q(a>) = 10 log q(u>). Assume, further, that the 
sample errors e, = (P, — Ri) are independent, zero-mean Gaussian 
variables, each having an rms value a me as. Since the Ri's are unknow- 
able, the e,'s cannot be observed directly; at best, their statistics can 
be inferred from an analysis of the observed errors, Ei = P, — Q,. 

Specifically, suppose that M is a sufficiently high polynomial order 
that, in the absence of measurement noise, all significant variations of 
Pi with i could be accommodated by the fitting function Q(u). In that 
case, all the Ej's would be due solely to measurement errors, the e/s. 
To a first approximation, we could then say the following: 8 The 2?,'s 
are zero-mean Gaussian variables, each having an ensemble mean- 
square value 

22 — M 

E'i = — -- — aLas; each i. (21) 

In addition, E^ is chi-square distributed, with 22 - M degrees of 
freedom and a mean given by eq. (21), i.e., ELs = E'f . Finally, the 
errors of major interest, e, = (Qi — R,), are zero-mean Gaussian 
variables, each having a mean square value 

€ fm of =tf-Ei 

= a Las! each i. (22) 

The pertinence of this theory to our data reductions is as follows: If 
the observed Ei's are essentially zero-mean and Gauss-distributed 



* Remember that the least -squares method is applied to E™ rather than to the 
quantity of interest, Emi*. The possibility thus exists of a small bias error in the results, 
which this step removes. 
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Fig. 3— Two examples of data records and their representations by zeroth-, second- 
and fourth-order polynomials. 



across the data ensemble, each with the same E~i ; and if E m s is chi- 
square distributed with 22 — M degrees of freedom and a mean of Ef , 
then the random-error model, as described above, can be assumed to 
apply. Accordingly, oiLas can then be estimated from measured error 
statistics using eq. (21), and the mean-square error between the fitted 
and true power gains can be estimated using eq. (22). 

The reduced data show that, when M = 0, the populations of E'^ 
and the E's do not have the properties of the random error model. 
This is to be expected, since a zeroth-order polynomial (horizontal 
line) cannot accommodate the true variations of faded power with 
frequency; hence, the E's are caused mostly by the inadequacy of the 
fitting function rather than by measurement noise. 
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model. 



For M = 2 and M = 4, however, the populations of i?Ls and the EJb 
do exhibit the properties of the random error model. For example, Fig. 
4 shows variations of Ei, E\ 2) and E\ 4) with i for M = 2 (see Table I). 
Each vertical bar spans the spread of values among the four data 
subsets, while the dashed lines give the values predicted by the random 
error model. 

Figure 5 gives even stronger evidence for the random error model 
when M > 2. The vertical axis represents Pr{£ rms > abscissa}, the 
vertical bars give the spread of this quantity among the four data 
subsets, and the dashed curves give the variations predicted by the 
random error model when a meas = 0.63 dB. For M = 0, the measured 
results compare very poorly with the predictions, for the reasons noted 
above, while for M = 2 and M = 4 they compare very well. Using eq. 
(22), the errors between Q, and Ri are found to have the rms values 
a e = 0.23 dB for M = 2 and o e = 0.29 dB for M = 4. 
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The rms modeling error, o e , is smaller for M = 2 because only three 
polynomial coefficients are derived from 23 recorded samples, leaving 
20 degrees of freedom for averaging out measurement noise. For M = 
4, five polynomial coefficients are derived, leaving 18 degrees of free- 
dom for noise averaging. We conclude that, for modeling purposes, M 
should be large enough to provide adequate curvature in q(w) (as 
evidenced by E,- and EJLu-populations that fit the random error 
model), but no larger (so as to minimize a e ). In view of this criterion, 
we have decided that M = 2 is the optimal polynomial order for the 
present data base. 

Table II gives computed ensemble averages for £L (the statistic 
minimized by least-squares methods) and E'^ for M = 0, 2, and 4. The 
ensemble in this cas e wa s the four merged data subsets. The close 
agreements between Eim* and £L for both values of M justify the 
approximation of Ei™ by EiLs in the least-squares derivation of the 
D m 's. 

The results in Table II also permit comparisons between the cases 
where the recorded data are adjusted, as described in Section 4.2.1, 
and not adjusted. The reduction in fitting error associated with these 
adjustments is small but distinct, indicating that the inferred calibra- 
tion differences are real. 

4.2.2.3 Coefficient Distributions. We have determined {Do, D\, Dz) for 
every record in the total set, using the least-squares fitting procedure 
described above. This effort has led to the set of cumulative probability 
distributions in Figs. 6, 7, and 8. In each case, the vertical bars span 
the results for the four data subsets, while the solid curve gives the 
result for the total set. In some cases, where the probability curve is 
steep and the spreads are small, the vertical bars can barely be 
distinguished. In all three figures, the results show good consistency 
among the four subsets and the total set. 

The distribution for Do, Fig. 6, is worth special note. From eq. (3), it 
is clear that Do can be interpreted as the power-gain ratio at the center 
of the channel (to = 0). Accordingly, it must always be positive, whereas 
D\ and D> can have either polarity. Also, since the power-gain statistics 
should be the same at each frequency, the decibel value of Do should 
have the same distribution as do the P,'s. For that reason, the abscissa 

Table II — Error statistics for 8640 records 





M = 


M 


= 2 


M 


= 4 


Measured 

Error 
Statistics 


Data No Adjust- 
Adjusted ment 


Data 
Adjusted 


No Adjust- 
ment 


Data 
Adjusted 


No Adjust- 
ment 


1^ in dB 2 
E^, in dB 2 


1.698 1.880 


0.361 
0.351 


0.486 
0.473 


0.314 
0.311 


0.373 
0.368 



POLYNOMIAL MODEL FOR MULTIPATH FADING 1213 



99.8 

99.5 

99 


Vs. 






98 








95 








90 


- 






Z 
UJ 

o 80 
cc 








CUMULATIVE PROBABILITY IN P 
O O O O 




| SPAN OF RESULTS FOR THE 
FOUR DATA SUBSETS 


-Do 












\ 


5 
2 


CUMULATIVE PROBABILITY FO 

P, (ALL i, TOTAL SET) 


R 


1 
0.5 


- 




^i 


0.1 




I 1 1 I 1 


1 1 



-41 



-37 



-33 



-25 -21 

10 log D OR Pj 



-5 



Fig. 6— Probability distributions for power-gain polynomial coefficient Do and re- 
corded, adjusted dB power-gain P,. 

in Fig. 6 is scaled in units of 10 log Do, and the distribution for all Pfa 
in the data base (24,820 records X 23 samples/record) is given as well 
(dashed curve). The agreement between the dashed and solid curves 
is excellent, as we should expect. 

It should also be noted that a straight line on the probability paper 
used here would imply a Gaussian (or normal) probability law. Thus, 
10 log Do is nearly normal, i.e., D is approximately log-normal. Since 
Aq = yfDo (eq. (7a)), this means that A is also approximately log- 
normal, as we have noted in Section III. 

4.2.3 Coefficients of the complex gain polynomial 

The method for computing the A„'s and £„'s from the D m 's for M 
= 2 was outlined in Section 2.2. The primary results are given by eqs. 
(6) and (7). The condition A > occurs over 58 percent of the data 
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records, in which cases | H c (w) | 2 can be matched to q(to) without 
requiring a second-order coefficient (i.e., A-i = B2 = 0). Over the 42 
percent of the records wherein A < 0, however, A2, £2 5* 0, and a small 
but nonzero fourth-order term, (A| + BDu 4 , exists in | H c (u) | 2 . Al- 
though the method of solution is designed to minimize this term, it is 
necessary to evaluate its effect on the polynomial fitting. 

We define the added distortion, 8, for a given record to be the rms 
decibel difference between | H c (u) \ 2 and q(io), where the rms averaging 
is over i. Thus, 



8 = \/ (10 log 



q(wi) + (A| + Bl)u 



q(ui) 



[ }' : 



(23) 



Clearly, 8 = for records in which A > 0. We have computed a 
cumulative distribution for nonzero 5 over 3600 records, specifically, 
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those records within the four merged subsets where A < 0. The result 
(Fig. 9) is that, for 99.4 percent of these records, 8 is below 0.23 dB, i.e., 
smaller than the rms error due to noise. Over all records, moreover, 8 
is less than 0.52 dB. 

Having established the accuracy of the procedure for computing the 
A n 's and B n 's, we have elected to ignore those nonzero values of A 2 
and B 2 that arise in 42 percent of the records. This step permits a 
simpler model, as noted previously, at some cost in accuracy. Justifi- 
cation for it is given in Section V. 

4.3 Deriving the coefficient statistics 

Assume now that the set (A , A h Bi) is derived, record by record, for 
the total data base. Using this large population of coefficient sets, it 
should be possible to obtain a mathematical description for the pdf 
p(A , Ai, Bi). The procedure used was outlined in Section 2.3, and the 
result was presented in Section III. In the following discussion, we 
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present the evidence for (i) the pdf of a„ [eqs. (14) to (16)]; (ii) the 
pdf's of A\ and B\ conditioned on a a ; {Hi) the functional variations of 
oa and ob with a a [eqs. (17) and (18)]; and (iv) the statistical indepen- 
dence of A] and Bi for given a . 

4.3.1 PDF for a 

The quantity a Q is a shifted, scaled version of 20 log An which, from 
eq. (7a), is identical to 10 log Do. The empirical distribution for this 
variate in Fig. 6 is represented in Fig. 10 by the circles. The dashed 
curve in Fig. 10 is based upon the pdf for a given by eqs. (15) and (16). 
We conclude, then, that our mathematical model for the distribution 
of 20 log Ao is an accurate one. 

4.3.2 Conditional PDF's for A^ and Si 

To derive the pdf's of A\ and B\, conditioned on a,„ the 24,820 data 
records were grouped into 11 subpopulations, each corresponding to a 
particular range of a„. For each of the central nine subpopulations, the 
a„-range has a width of one-half a standard deviation, thereby balanc- 
ing the objectives of good resolution (narrow a„-range per subpopula- 
tion) and sample size sufficiency (many records per subpopulation). 

Within each subpopulation, we computed the mean, standard devia- 
tion and cumulative distribution for Ai, and similarly for Bi. Since the 
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polarity of Ai was chosen randomly in each data record (a reaction to 
the absence of phase response data, as noted in Section 2.2), the mean 
for A\ was very close to zero in every subpopulation. 

We also randomized the polarity of every reduced value of B\ before 
analyzing the B\ statistics in each subpopulation. The motive for this 
polarity randomization is that it simplifies the statistical modeling. 
Before randomization, the conditional mean of B\ was not close to zero 
in every subpopulation, nor was its conditional distribution Gauss-like; 
after randomization, B\ became a zero-mean Gauss-like variable in 
every subpopulation. Although this simplification entailed some doc- 
toring of the derived Bi's, it should be noted that the polarity of B Y is 
immaterial in any practical system context. For this quantity repre- 
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sents the slope, at band center, of the real part of H c (<o), and whether 
the slope is upward to the right (Bi negative) or upward to the left (jBj 
positive) has no effect on either detection behavior or the realizability 
of equalizer circuits. The same cannot be said for the polarity of A\ t 
which can influence the stability of some equalizer designs. 

Table III gives some pertinent data for the 11 subpopulations 
analyzed. For each one, identified by a number and an a Q -range, the 
table gives the number of records contained therein and the standard 
deviations for A\ and B\. Because of the polarity randomizations, the 
means for A\ and B x are essentially zero in every case. 

Some distributions for A\/oa and B\/ob are indicated in Figs. 11 and 
12, respectively. Data are given for only five of the subpopulations, but 
all of the nine central subpopulations exhibit the Gauss-like distribu- 
tions evident here. That is, A\/oa and Bi/ob act, in every case, like 
zero-mean, unit-variance Gaussian variables. We can therefore invoke 
the following models: 

p A (A! I a ) = —= exp j - - 2 ' \ (24) 



and 



1 f 1 B 2 

p B (B, I a n ) = -j= exp] - - -2-^- \ . (25) 



4.3.3 a A (a ) and a B (a ) 

Let us associate each oa in Table III with the midvalue of the ar- 
range for which it is computed and do the same for each ob- (This 
procedure has no meaning, of course, for the semi-infinite ranges of 
the first and last subpopulations.) The result is that the circles and 
crosses in Fig. 13 can be regarded as samples of the functions oa(o,o) 
and ob{cIo), respectively. The solid curves provide accurate fits to these 
data that are described mathematically by eqs. (17) and (18). 

Table III — Data for the 1 1 subpopulations 



Subpop. 




Population 


o A 


OB 


Number 


a„-Interval 


Size 


(ns) 


(ns) 


1 


< -2.25 


515 


0.154 


0.126 


2 


[-2.25, 1.75) 


660 


0.153 


0.125 


;i 


[-1.75, -1.25) 


1423 


0.141 


0.121 


4 


[-1.25, -0.75) 


2878 


0.176 


0.140 


5 


[-0.75, -0.25) 


4147 


0.240 


0.186 


6 


[-0.25, 0.25] 


4925 


0.278 


0.175 


7 


(0.25, 0.75] 


4340 


0.370 


0.185 


8 


(0.75, 1.25] 


3522 


0.440 


0.239 


9 


(1.25, 1.75] 


1933 


0.503 


0.233 


10 


(1.75, 2.25] 


371 


0.576 


0.218 


11 


>2.25 


J106 
Total: 24,820 


0.656 


0.158 
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Fig. 11 — Probability distributions for A 1/0,4, conditioned on a„, for five subpopula- 

tions. 



4.3.4 Statistical independence of A ^ andB^ 

Formal statistical testing has shown that A\ and Si within each 
subpopulation are not mutually independent. We have assumed such 
independence anyway because of the resulting simplification of the 
model. The validity of this approach rests upon how it influences the 
outcome of system studies. The next section takes up this issue and 
gives justification for treating A\ and B\, conditioned as a , as inde- 
pendent variables. 

V. ACCURACY OF THE MODEL 

The multipath fading model presented here consists of (i) a poly- 
nomial representation for H c (u); (ii) a mathematical formula for the 
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joint coefficient pdf; and [Hi) a formula for Tm, the number of fading 
seconds per heavy-fading month. We now discuss the accuracy of each 
of these constituents. 

5. 1 The polynomial representation 

The first step in using the polynomial representation consists of 
fitting q(io), eq. (3), to the recorded power gain data. This was done, 
with M = 2, for each of 24,820 records. The error measured in each 
record is the rms decibel difference (.E rms ) between the recorded and 
fitted power gains. Although the error between q(u) and the true 
underlying response cannot be estimated from isolated records, its 
statistics can be inferred by analyzing the population of Enos over 
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thousands of records. Such an analysis* had led to the following 
conclusion: The decibel difference between the fitted and true power 
gains is approximately Gauss-distributed at each frequency, is due 
almost entirely to measurement noise (rather than inadequate curva- 
ture in the form of q(u)), and has an rms value of 0.23 dB. This 
corresponds to an rms error of 2.7 percent in the gain magnitude, \H C \, 
which is sufficiently low for purposes of modeling. 

The next step in the polynomial representation is to derive the A- 
and ^-coefficients in eq. (2) from the D-coefficients in eq. (3). For 58 
percent of the data records, coefficient sets (Ao, A\, B\) exist which 
provide a perfect match between | H c (u>) | 2 and q(co). The one flaw in 
these cases is the unresolvable ambiguity (resulting from noncoherent 
measurements) in the polarity ofAi. 

For the remaining 42 percent of the records, the method for choosing 
A's and B's leads, in addition, to nonzero A 2 and B 2 and a resulting 
fourth-order term in \H c (<o) | 2 not present in q(w). The differences 
between | H c {u>) \ 2 and q(u>) for these records were analyzed, and found 
to be minor (see Section 4.2.3 for more details). 

Finally, all nonzero values computed for A2 and B2 have been 
discarded to simplify the model. To evaluate the validity of this step, 
we performed an outage analysis for several different modulations, 
specifically, 4-psk, 8-psk, 16-psk and 16-qam. For each modulation, 
cosine rolloff spectral shaping was assumed, and the rolloff factor (a) 
and channel bandwidth (W) were treated as variables. In every case, 



* It is noteworthy that the distribution for £L derived here using a second-order 
q(u) is nearly identical to that obtained for the three-path model of Rummler (Ref. 5). 
The slightly lower mean obtained here is due solely to the data adjustments described 

in Section IV. 
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we specified the symbol rate to be W/(l + a), thereby matching the 
truncation bandwidth of the modulation to the channel bandwidth. 

The outage analysis performed for each modulation used W = 20, 
30, and 40 MHz and a-values between 0.1 and 1.0. For each of the 
resulting cases, the quantity "outage seconds/heavy-fading month" 
was estimated using the following procedure: Regions of the coefficient 
set (Ao, Ai, B\, A 2 , B 2 ) for which the detected data eye is closed were 
derived analytically;* the number of data records for which {A , - - • 
B 2 ) lies within these so-called outage regions was counted; and the 
fraction of these records was multiplied by T M in eq. (12). [We assumed 
c = 1 (average terrain), F = 6 (GHz), and d = 26.4 (miles), resulting in 
a value for T M of 12,144 s]. 

This computation was done twice for each combination of modula- 
tion, W and a. In one computation, nonzero values of A 2 and B 2 were 
acknowledged in comparing each (Ao, • • • B 2 ) to the outage regions; in 
the other, A 2 and B 2 in each record were taken to be zero. 

The numerical results span a range from 5 to 8000 seconds of outage. 
The results obtained by acknowledging nonzero (A 2 , B 2 ) pairs and the 
results obtained by ignoring them differ by less than 10 percent for all 
cases treated. The simplification of omitting A 2 and B 2 from the model 
is thus found to be acceptable for purposes of making outage predic- 
tions. 

5.2 The coefficient PDF 

In deriving the mathematical result for p(a„, A\, B\), Section III, 
numerous approximations and assumptions were made, e.g., indepen- 
dence between A\ and B\ for given a a ; Gaussian conditional pdf 's for 
A] and B U ' and the functions for p a (a ), o A (a ), and ob(cio). The ultimate 
test of the result is whether it leads to the same performance predic- 
tions as those based upon the actual population of coefficient sets. 
Accordingly, outage seconds for the combinations of modulation, W, 
and a cited above were obtained using a separate method of calcula- 
tion: In addition to counting data records over the coefficient outage 
regions, the model distribution p(a a , A\, B\) was integrated over these 
regions. 

The main findings are as follows: When the outage calculated by 
counting data records is greater than 70 seconds per heavy-fading 
month, the outage calculated by integrating p(a a , A h Bi) differs by 16 
percent or less; for outages less than 70 seconds, the difference is 50 
percent or less. This is an excellent level of agreement for practical 
purposes. 



* The closed-eye criterion for "outage" leads to conservative performance measures 
but has the virtue of not depending on system fade margin or specified bit error rate; it 
also simplifies the outage analysis considerably. 
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The above findings vindicate the modeling assumption that Ai and 
B\ are independent for given a„. It also confirms the accuracy of the 
functions used for p a (a ), oa(cIo), and ob(cIo). It is noteworthy that, 
when outage is calculated using a purely Gaussian function for p a (a ) 
[i.e., z(a ) = in eq. (16)], the agreement with counted data records is 
poor in several cases. In short, the refinement in p a (a ) used here is 
important to the accurate use of the model in system outage predic- 
tions. 

5.3 Expression for T» 

The expression given for Tm in Section III assumes only that this 
quantity is in a fixed proportion to the multipath occurrence factor of 
Barnett. 7 By computing the fading time and occurrence factor for the 
present data base, we have obtained the empirical results of eq. (12). 
This result is unverified, although recent studies have indicated that 
it is reliable to within a factor or two. 9 

There is further evidence supporting eq. (12) for Tm, as well as eq. 
(1) for p a (a )' Using these two formulas, plus (14) and (16), we can 
show that deep fades observe the probability law 

Pr{A < L) = 1.7 rL 2 \ (-10 log L 2 > 26 dB). (26) 

This law applies to heavy-fading months, and r is the multipath 
occurrence factor used to derive Tm- This result differs from that of 
Barnett 7 by the factor 1.7, which represents very good agreement. 

VI. DISCUSSION AND CONCLUSION 

We have presented a statistical model for describing the microwave 
channel distortions classified as multipath fading. The model is defined 
by the functions of eqs. (11) through (18) and by the 12 numerical 
constants contained therein. We do not know whether the functional 
structure is valid for other paths, microwave frequencies, seasons, etc., 
and, if it is, whether or how the numerical constants vary with these 
conditions. The way to resolve these questions is through more mul- 
tipath experiments. 

A major virtue of the polynomial model suggested by eq. (2) is the 
simplicity it affords in analyzing the response of a fading channel to 
arbitrary modulations. This is largely because the required polynomial 
order for H c (io) in channelized microwave systems is generally quite 
low. Ample support for this conclusion is given by the data reductions 
and system analysis results reported here. 

Another useful feature of the polynomial model is that it suggests 
the form of an adaptive equalizer response, i.e., 1/H c (u), which is 
easy to realize when the complex zeros of H c (co) have negative real 
parts. With this in mind, the polarity of Ai becomes a very important 
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issue. For when A\ > 0, the realizable adaptive response [A a — wBi + 
ja \Ai | ] _1 is highly effective against multipath fading. 10 When A\ < 0, 
however, such a response eliminates amplitude distortion at the cost 
of increased delay distortion. In lieu of explicative data, the present 
model assumes that A\ < precisely half the time. Whether this is 
true or not has important bearing on the achievement of successful 
equalizer approaches. Coherent multipath measurements would help 
to resolve the existing polarity uncertainty in the model. 

The new model provides a useful starting point for the design of 
subsequent experiments. What it reveals about the depth, shape, and 
statistics of multipath fading responses can facilitate design choices 
for the channel probing signal, receiver processing, data acquisition 
and reduction strategies, and other experiment features. 
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